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Abstract 

For a regular chain R, we propose an algorithm which computes the (non-trivial) 
limit points of the quasi-component of R, that is, the set W{R) \ W{R). Our proce- 
dure relies on Puiseux series expansions and does not require to compute a system of 
generators of the saturated ideal of R. We focus on the case where this saturated ideal 
has dimension one and we discuss extensions of this work in higher dimensions. We 
provide experimental results illustrating the benefits of our algorithms. 

1 Introduction 

The theory of regular chains, since its introduction by J.F. Ritt PH, has been applied 
successfully in many areas including parametric algebraic systems [TU], differential sys- 
tems [m 131 [17] , difference systems [16] , intersection multiplicity [21] , unmixed decomposi- 
tions [18j and primary decomposition |27j of polynomial ideals, cylindrical algebraic decom- 
position [5], parametric [3T] and non- parametric [S] semi-algebraic systems. Today, regular 
chains are at the core of algorithms for triangular decomposition of polynomial systems, 
which are available in several software packages [THl (2^1 130] ■ Moreover, these algorithms 
provide back-engines for computer algebra system front-end solvers, such as Maple's solve 
command. 

One of the algorithmic strengths of the theory of regular chains is its regularity test 
procedure. Given a polynomial p and a regular chain R, both in a multivariate polynomial 
ring k[Xi, . . . , Xn] over a field k, this procedure computes regular chains . . . , i?e such 
that Ri, . . . ,Re is a decomposition of R in some technical sense and for each 1 < i < e 
the polynomial p is either null or regular modulo the saturated ideal of Ri. Thanks to the 
D5 Principle |13j . this regularity test avoids factorization into irreducible polynomials and 
involves only polynomial GCD and resultant computations. 

One of the technical difficulties of this theory, however, is the fact that regular chains 
do not fit well in the "usual algebraic-geometric dictionary" (Chapter 4, [12J. Indeed, 
the "good" zero set encoded by a regular chain i? is a constructible set W{R), called the 
quasi-component of R, which does not correspond exactly to the "good" ideal encoded by 
R, namely sat(i?), the saturated ideal of R. In fact, the affine variety defined by sat(i?) 
equals W{R), that is, the Zariski closure of W{R). 

For this reason, a decomposition algorithm, such as the one of M. Kalkbrener [TS] 
(which, for an input polynomial ideal I computes regular chains Ri, . . . , Re such that 
equals the intersection of the radicals of the saturated ideals of . . . , Re) can not be seen 
as a decomposition algorithm for the variety V{I). Indeed, the output of Kalkbrener's 
algorithm yields V{I) = W{Ri) U • • • U W{Re) while a decomposition of the form V{I) = 
W{Ri) U • • • U W{Rf) would be more exphcit. 

^The radical of the saturated ideal of R is equal to the intersection of the radicals of the saturated ideals 
of Rl,...,Re. 
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Kalkbrener's decompositions, and in fact all decompositions of differential ideals [HI El 
I17j raise another notorious issue: the Ritt problem, stated as follows. Given two regular 
chains (algebraic or differential) R and S, check whether the inclusion of saturated ideals 
sat(i?) C sat(S') holds or not. In the algebraic case, this inclusion can be tested by com- 
puting a set of generators of sat(i?) , using Grobner bases. In practice, this solution is too 
expensive for the purpose of removing redundant components in Kalkbrener's decomposi- 
tions and only some criteria are applied [20]. In the differential case , there has not even 
an algorithmic solution. 

In the algebraic case, both issues would be resolved if one would have a practically 
efficient procedure with the following specification: for the regular chain R compute regular 
chains Ri, . . . ,Re such that we have W{R) = W{Ri) U • • • U W{Re). If in addition, such 
procedure does not require a system of generators of sat(i?), this might suggest a solution 
in the differential version of the Ritt problem. 

In this paper, we propose a solution to this algorithmic quest, in the algebraic case. 
To be precise, our procedure computes the non-trivial limit points of the quasi-component 
W{R), that is, the set lim(W^(i?)) := W{R)\ W{R). This turns out to be W{R) n V{hR), 
where V{hR) is the hypersurface defined by the product of the initials of R. We focus on 
the case where the saturated ideal of R has dimension one. In Section [lO) we sketch a 
solution in higher dimension. 

When the regular chain R consists of a single polynomial r, primitive w.r.t. its main 
variable, one can easily check that \im{W{R)) = V{r, hn) holds. Unfortunately, there is no 
generalization of this result when R consists of several polynomials, unless R enjoys remark- 
able properties, such as being a primitive regular chain |20j . To overcome this difficulty, 
it becomes necessary to view i? as a "parametric representation" of the quasi-component 
W{R). In this setting, the points of \im{W{R)) can be computed as limits (in the usual 
sense of the Euclidean topology of sequences of points along "branches" (in the sense of 
the theory of algebraic curves) of W{R) . It turns out that these limits can be obtained 
as constant terms of convergent Puiseux series defining the "branches" of W{R) in the 
neighborhood of the points of interest. 

Here comes the main technical difficulty of this approach. When computing a particular 
point of lim(W^(i?)), one needs to follow one branch per defining equation of R. Following a 
branch means computing a truncated Puiseux expansion about a point. Since the equation 
of R defining a given variable, say Xj, depends on the equations of R defining the variables 
Xj^i, Xj^2, ■ ■ ■, the truncated Puiseux expansion for Xj is defined by an equation whose 
coefficients involve the truncated Puiseux expansions for Xj_i, Xj_2, .... 

From Sections [4] to [Sj we show that this principle indeed computes the desired limit 
points. In particular, we introduce the notion of a system of Puiseux parametrizations 
of a regular chain, see Section [4j This allows to state in Theorem [3] a concise formula 
for lim(VF(i?)) in terms of this latter notion. Then, we estimate to which accuracy one 
needs to effectively compute such a system of Puiseux parametrizations in order to deduce 
\im{W{R)), see Theorem [6] in Section [t] 

In Section |9j we report on a preliminary implementation of the algorithms presented 
in this paper. We evaluate our code by applying it to the question of removing redundant 
components in Kalkbrener's decompositions and observe the benefits of this strategy. 

In order to facilitate the presentation of those technical materials, we dedicate Section [3] 
to the case of regular chains in 3 variables. Section[2]briefly reviews notions from the theories 

^This identification of tfie closures of W{R) in Zariski topology and the Euclidean topology holds when 
k is C. 
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of regular chains and algebraic curves. We conclude this introduction with a detailed 
example. 

Consider the regular chain R = {ri,r2} C k[Xi,X2,X3] with ri = XiX|+X2 + l,r2 = 
{Xi + 2)XiX^ + {X2 + 1)(X3 + 1). Then, we have Hr = Xi{Xi +2). To determine lim(i?), 
we need to compute Puiscux scries expansions of ri about Xi = and Xi = —2. We start 
with Xi = 0. The two Puiseux expansions of ri about Xi = are: 

[X,=T,X2 = ^^^ + 0{T^)], 
[X,=T,X2 = =^±^+0{T^)]. 

The second expansion does not result in a new limit point. After, substituting the first 
expansion into r2, wc have: 

r'2 = r2{Xi =T,X2^ + 0(T2), X3) 

= T {{T + 2) Xs^ + {0{T^)-5T^ - 2 T - 1 + l) (X3 + 1)) . 

Now, we compute Puiseux series expansions of r2 which are 

[T = T,X3 = 1 + T + CI(T2)], 
[T = T,X3 = -1/2- l/4r + 0(T2)]. 

So the regular chains {Xi , X2 + 1, X3 — 1} and {Xi , X2 + 1, X3 + 1 /2} give the limit points 
of W{R) about Xi = 0. 

Next, we consider Xi = —2. We compute Puiseux series expansions of ri about the 
point Xi = —2. We have: 

[Xi=T-2,X2 = l + l/3T + 0(r2)], 
[Xi = T-2,X2 = -1/2- 1/12 T + 0(r2)]. 

After substitution into r2, we obtain: 

r'i2 =r2(Xi =r-2,X2 = 1 + 1/3 r + 0(r2),X3) 
= (T - 2) TXs^ + (2 + 1/3 T + 0{T'^)) {X3 + 1) 
r'22 =r2{[Xi=T-2,X2 = -1/2- 1/12 r + 0(T2)) 
= (T - 2) TXs^ + (1/2 - 1/12 T + 0(T2)) {X3 + 1) . 

So those Puiseux expansions of r'12 and r22 about T = which result in a limit point are 
as follows: 

i) for r',2: [T = T,Xs = ^+ 0{T^)] 

ii) for r'22: [T = T,Xs = + 0(T2)] 

Thus, the limit points of R about the point Xi = —2 can be represented by the regular 
chains {Xi + 2, X2 - 1, X3 + 1} and {Xi + 2, X2 + 1/2, X3 + 1}. 

One can check that a triangular decomposition of the system R U {Xi} is {X2 + 1, Xi} 
and, thus, does not yield \im{W{R)) n V{Xi), but in fact a superset of it. 



3 



2 Preliminaries 



This section is a brief review of various notions from the theories of regular chains, alge- 
braic curves and topology. For these latter subjects, our references are the textbooks of 
R.J. Walker [3S], G. Fischer [IS] and J. R. Munkres [53]. The notations and hypotheses 
introduced in this section are used throughout the sequel of the paper. 

Multivariate polynomials. Let k be a field which is algebraically closed. Let Xi < 
■ ■ ■ < Xs he s > 1 ordered variables. We denote by k[Xi, . . . , Xs] the ring of polynomials 
in the variables Xi,...,Xs and with coefficients in k. For a non-constant polynomial 
p € k[Xi, . . . , Xs], the greatest variable inp is called main variable of p, denoted by mvar(p), 
and the leading coefficient of p w.r.t. mvar(p) is called initial of p, denoted by init(p). 

Zariski topology. We denote by the ajfine s-space over k. An ajfine variety of A'* is 
the set of common zeroes of a collection F C k[Xi, . . . ,Xs] of polynomials. The Zariski 
topology on A" is the topology whose closed sets are the affine varieties of A* . The Zariski 
closure of a subset V7 C A* is the intersection of all affine varieties containing W. This is 
also the set of common zeroes of the polynomials in k[Xi, . . . , X^] vanishing at any point 
of T^. 

Relation between Zariski topology and the Euclidean topology. When k = C, the affine 
space A* is endowed with both Zariski topology and the Euclidean topology. The basic open 
sets of the Euclidean topology are the balls while the basic open sets of Zariski topology 
are the complements of hypersurfaces. A Zariski closed (resp. open) set is closed (resp. 
open) in the Euclidean topology on A*. The following properties emphasize the fact that 
Zariski topology is coarser than the Euclidean topology: every nonempty Euclidean open 
set is Zariski dense and every nonempty Zariski open set is dense in the Euclidean topology 
on A'*. However, the closures of a constructible set in Zariski topology and the Euclidean 
topology are equal. More formally, we have the following (Corollary 1 in LIO of |23j) key 
result. Let C A" be an irreducible affine variety and U C V he open in the Zariski 
topology induced on V. Then the closure of U in Zariski topology and the closure of U in 
the Euclidean topology are both equal to V. 

Limit points. Let {X, r) be a topological space. A point p £ X is a limit of a sequence 
{xn,n G N) of points of X if, for every neighborhood U of p, there exists an N such that, 
for every n > N, we have Xn G U; when this holds we write lim„_j.oo Xn = P- If X is a 
Hausdorff space then limits of sequences are unique, when they exist. Let S C X he a 
subset. A point p g X is a limit point of S if every neighborhood of p contains at least one 
point of S different from p itself. Equivalently, p is a limit point of S if it is in the closure 
of 5 \ {p}. In addition, the closure of S is equal to the union of S and the set of its limit 
points. If the space X is sequential, and in particular if X is a metric space, the point p is 
a limit point of S if and only if there exists a sequence {Xn, n £ N) of points of 5 \ {p} with 
p as limit. In practice, the "interesting" limit points of S are those which do not belong to 
S. For this reason, we call such limit points non-trivial and we denote by lim(5) the set of 
non-trivial limit points of S. 

Regular chain. A set R of non-constant polynomials in k[A'i, . . . , Xg] is called a triangular 
set, if for all p,q £ R with p ^ q we have mvar(p) ^ mvar(g). For a nonempty triangular set 
R, we define the saturated ideal sat(i?) of R to be the ideal (i?) : h^, where Hr is the product 
of the initials of the polynomials in R. The empty set is also regarded as a triangular set, 
whose saturated ideal is the trivial ideal (0). From now on, R denotes a triangular set of 
k[Ari, . . . ,Xs\. The ideal sat(i?) has several properties, in particular it is unmixed [4]. We 
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denote its height by e, thus sat(i?) has dimension s—e. Without loss of generality, we assume 
that k[Xi, . . . , Xs_f,] n sat(i?) is the trivial ideal (0). For all 1 < i < e, we denote by the 
polynomial of R whose main variable is Xi^^-e and by hi the initial of r^. Thus hn is the 
product hi ■ ■ ■ he- We say that i? is a regular chain whenever R is empty or {ri, . . . , r^-i} 
is a regular chain and he is regular modulo the saturated ideal sat({ri, . . . ,re-i}). The 
regular chain R is said strongly normalized whenever hn G k[Xi, . . . ,Xs_e] holds. If R is 
not strongly normalized, one can compute a regular chain N which is strongly normalized 
and such that sat(i?) — sat(A^) and V{hN) — V{hij) both hold, where hji is the iterated 
resultant of hji w.r.t R. See 7 . 

Limit points of the quasi-component of a regular chain. We denote by W{R) :— V{R) \ 
V{hii) the quasi- component of R, that is, the common zeros of R that do not cancel h^^. 
The above discussion implies that the closure of W{R) in Zariski topology and the closure 
of W{R) in the Euclidean topology are both equal to V {sat{R)) , that is, the affine variety 
of sat(i?). We denote by W{R) this common closure. We call limit points of W{R) the 
elements of hm(W^(i?)). 

Rings of formal power series. Recall that k is an algebraically closed field. From now 
on, we further assume that k is topologically complete. Hence k may be the field C of 
complex numbers but not the algebraic closure of the field Q of rational numbers. We 
denote by k[[Xi, . . . , X^]] and \i{Xi, . . . , Xg) the rings of formal and convergent power 
series in Xi, . . . , Xg with coefficients in k. Note that the ring k(Xi, . . . , Xg) is a subring 
of k[[Xi, . . .,Xs]]. When s = 1, we write T instead of Xi. Thus k[[T]] and k(T) are 
the rings of formal and convergent univariate power series in T and coefhcients in k. For 
/ G k[[Xi, . . . , Xg]], its order is defined by 

r min{d I /(,) ^ 0} if /^O, 
\ oo if./ = 0. 

where /(^^ is the homogeneous part of / in degree d. Recall that k[[Xi, . . . , Xg]] is topologi- 
cally complete for Krull Topology and that k(Xi, . . . , Xg) is a Banach Algebra for the norm 
defined by || / = Se |ae|/o'' where / = l^eaeX" G k[[Xi, . . .,Xg]] and p = (pi, . . . , pg) e 
Mf,Q. We denote by Mg the only maximal ideal of k[[Xi, . . . , Xg]], that is, 

= {/ek[[Xi,...,X,]] |ord(/)>l}. 

Let / e k[[Xi, . . . , Xg]] with / ^ 0. Let fc e N. We say that / is (1) general in Xg if / 7^ 
mod Mg-i, (2) general in Xg of order k if we have ord(/ mod A^s-i) = k- 

Formal Puiseux series. We denote by k[[T*]] — [J^^i k[[r"]] the ring of formal Puiseux 
series. For a fixed ip € k[[T*]], there is an n e N>o such that p E k[[T"]]. Hence 
p = X]m=o '^rn.T^ , wlicrc G k. We call order of p the rational number defined by 
ord(v3) = min{^ | ^ 0} > 0. We denote by k((r*)) the quotient field of k[[r*]]. 

Convergent Puiseux series. Let p e C[[T*]] and n e N such that p = /(T") with 
/ G C[[r]] holds. We say that the Puiseux series p is convergent if we have / S C{T). 
Convergent Puiseux series form an integral domain denoted by C{T*); its quotient field is 
denoted by C{{T*)). For every p> e C((r*)), there exist n G Z, r e N>o and a sequence of 
complex numbers a„, On+i, ari+2, ■ • ■ such that we have 

CO 

p — \^ OmT^ and a„ ^ 0. 
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Then, we define ord((p) = 

Puiseux Theorem. If k has characteristic zero, the field k((r*)) is the algebraic closure 
of the field of formal Laurent series over k. Moreover, if k = C, the field C((T*)) is 
algebraically closed as well. From now on, we assume k = C. 

Puiseux expansion. Let B = C{{X*)) or C{{X*)). Let / G M[Y], where d := deg(/, Y) > 0. 

Let h := lc(/, y). According to Puiseux Theorem, there exists (pi E M, i = 1, . . . ,d, such 
that ^ = (Y — ifi) ■ ■ ■ (Y — ifd). We call ipi, . . . ,ip(i the Puiseux expansions of / at the 
origin. 

Puiseux parametrization. Let / € C{X) [Y]. A parametrization of / is a pair (iplT), ip{T)) 
of elements of C(T) for some new variable T, such that (1) f{iJj{T).ip{T)) = holds in C{T), 
(2) wc have < oTd{ijj{T)), and (3) ip{T) and (p{T) are not botli in C. The parametrization 
{■iIj{T), f{T)) is irreducible if there is no integer fc > 1 such that both il){T) and (^(T) are in 
C(r'^). We call an irreducible parametrization {ip{T), if{T)) of / a Puiseux parametrization 
of /, if there exists a positive integer ^ such that ip{T) — T'' . The index <r is called the 
ramification index of the parametrization (T*^, Lp{T)). It is intrinsic to / and q < deg(/, Y). 
Let Zi,...,Zt; denote the primitive roots of unity of order <; in C. Then ip{ziX^/''), for 
z = 1, . . . , are ? Puiseux expansions of /. 

We conclude this section by a few lemmas which are immediate consequences of the 
above review. 

Lemma 1. We have: \im{W{R)) = W{R) fl V{hii). In particular, \im{W{R)) is either 
empty or an affine variety of dimension s — e — 1 . 

Lemma 2. If R is a primitive regular chain, that is, if R is a system of generators of its 
saturated ideal, then we have \iiii{W{R)) — V{R) Ci V{hii). 

Lemma 3. If N is a strongly normalized regular chain such that sat(i?) = sat(A^) and 
V{hN) = V{hR) both hold, then we have \im(W{R}) C \im{W{N)). 

Lemma 4. Let x £ such that x ^ W{R). Then x G \im{W{R)) holds if and only if 
there exists a sequence {an, n e N) of points in A* such that an € W{R) for alln gN and 
lim^_^QQ an — X. 

Lemma 5. Recall that R writes {ri, . . . , re}. If e > 1 holds, writing R' = {ri, . . . , rg-i} 
and r = re, we have 

lim{W{R'Ur)) C \im{W{R')) D lim{W{r)). 

Lemma 6. Let e C((r*)) and let p/q G Q be the order of if. Let (Q!„,n G N) be a 

sequence of complex numbers converging to zero and let N be a positive integer such that 
{ip{an),n > N) is well defined. Then, if p/q < holds, the sequence ((p(a„),n > N) escapes 
to infinity while if p/q > 0, the sequence (<^(a„),n > N) converges to the complex number 

3 Basic techniques 

This section is an overview of the basic techniques of this paper. This presentation is meant 
to help the non-expert reader understand our objectives and solutions. In particular, the 
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results of this section are stated for regular chains in three variables, while the statements 
of Sections |4] to [8] do not have this restriction. 

Recall that i? C C [Xi , . . . , X^] is a regular chain whose saturated ideal has height 
1 < e < s. As mentioned in the introduction, we mainly focus on the case e = s — 1, that 
is, sat(-R) has dimension one. 

Lemma [l] and the assumption e = s — 1 imply that \im{W{R)) consists of finitely many 
points. 

We further assume that R is strongly normalized, thus we have hji lies in C[Xi]. 

Lemma[2]and the assumption hn G C[Xi] imply that computing \im{W{R)) reduces to 
check, for each root a G C of ha whether or not there is a point x e lim(W^(i?)) whose 
Xi-coordinate is a. Without loss of generality, it is enough to develop our results for the 
case a = 0. Indeed, a change of coordinates can be used to reduce to this latter assumption. 

We start by considering the case n — 2. Thus, our regular chain R consists of a 
single polynomial ri £ C[Xi,X2] whose initial hi satisfies /ii(0) — 0. Lemma [t] provides 
a necessary and sufficient condition for a point of {a, P) £ A^, with a = 0, to satisfy 
(a,/3) elim(T^({ri})). 

Let d be the degree of ri in X2 ■ Applying Puiseux Theorem, we consider ipi, . . . ^ ipd G 
C{{Xi)) such that the following holds 

= (X2 - ^1) • • • (X2 - ip,i) (1) 

hi 

in C{{Xi))[X2]. We assume that the series ipi, . . . ,ipd are numbered in such a way that 
each of ipi, . . . ,ipc has a non-negative order while each of (pc+i, ■ ■ ■ has a negative order, 
for some c such that < c < d. 

Lemma 7. With hi{0) = 0, for all /3 € C, the following two conditions are equivalent 
(i) (0,/3) e lim(VF(ri)) holds, 

(ii) there exists 1 < J < c and a sequence (a„,n G N) of complex numbers such that the 
sequence {ipj{an), n G N) is well defined, we have hi{an) ^ for all n £ N and we we 
have 

lim an — and lim (/3,(q;„) — j3. 

n— >C30 n^oo 

Proof. We first prove the implication (ii) ^ (i). Equation ([T]) together with (ii) implies 
{an,(pj{an)) G V{ri) for all n G N. Since we also have (a„, (y9j(Q!„)) ^ V{hi) for all n G N 
and lim„_j.oo (ctn, fjic^n)) = (0, /?), we deduce (?), thanks to Lemma|4] 

We now prove the implication (i) (ii). By Lemma |4j there exists a sequence 
((Q!„,/3„),n G N) in such that for all n G N we have: (1) 7^ 0, (2) ri(Q!„,/3„) = 0, 

and (3) lim„_^oo (^n,/?™) = Since lim„^oo = 0, each series (pi(a„), . . . , V'd(an) is 

well defined for n larger than some positive integer N. Hypotheses (1) and (2), together 
with Equation Q, imply that for all n > iV the product 

(/3„ - ipiia„)) ■■■{fin - Vc{an))il3n - ipc+lictn)) ' ' ' - ipd{an)) 

is 0. Since lim„_5.oo (3n — fi, and by definition of the integer c, each of the sequences (/3„ — 
ipi{an)), . . . , iPn-^c{an)) couverges while each of the sequences (/3„-(pc+i(an)), ■ • • , (/3n- 
^d{ctn)) escapes to infinity. Thus, for n large enough the product (/3„ — Lpi{an)) ■ ■ ■ — 
^cic^n)) is zero. Therefore, one of sequences — <^i(q!„)), . . . , (/?„ — Lpc{<^n)) converges to 
and the conclusion follows. □ 
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Lemmas [6] and [T] immediately imply the following. 
Proposition 1. With ft-i(O) = 0, for all 13 E C, we have 

(0,/3)elim(M/(ri)) ^ /3 e {^i(0), . . . , (^^(0)}. 

Next, we consider the case n = 3. Hence, our regular chain R consists of two polynomials 
ri e C[Xi,X2] and r2 € C[Xi, X2, X^] with respective initials hi and /i2- We assume that 
is a root of the product /iift.2 and we are looking for all /3 € C and all 7 e C such that 
(0,/3,7)elim(W^(ri,r2)). 

Lemma [5] tells us that (0,;9,7) g lim(M^(ri, r2)) implies (0, /3) £ lim(iy(ri)). This 
observation together with Proposition [T] yields immediately the following. 

Proposition 2. With hi{0) = and /i2(0) 7^ 0, assuming that ri is primitive over C[Xi], 
for all f3 Cz C and all 7 G C, we have 

(0,/3,7)elim(W^(ri,r2)) ^ (0, /3, 7) € F(ri, r2). 

We turn now our attention to the case /ii(0) = /i2(0) = 0. Since (0,/3) G lim(M^(ri)) 
is a necessary condition for (0,/3,7) € \im{W{ri,r2)) to hold we apply Proposition [T] and 
assume (3 S {(^i(0), . . . , </?c(0)}. Without loss of generality, we further assume /3 = 0. 
For each 1 < j < c, such that ^Pj{0) = holds, we define the univariate polynomial 
/|eC((Xi*))[^3]by 

/|(Xi,X3) -r2(Xi,^,(Xi),X3). (2) 

Let b be the degree of . Applying again Puiseux theorem, we consider ipi, . . . ,ipb S 
C{{X'()) such that the following holds 

^ - (X3 - Vi) • • • (^3 - ^b) (3) 

'12 

in C((X^))[X3]. We assume that the series i/^i, . ..jtpb are numbered in such a way that 
each of V'l , ■ • ■ , V'a has a non-negative order while each of ipa+i, ■ • ■ i V'fc has a negative order, 
for some a such that < a < 6. 

Lemma 8. For all 7 G C, the following two conditions are equivalent, 
(i) (0,0,7) e lim(M/(ri,r2)) 

(m) </iere exist integers j, fc mt/j 1 < j < c anrf 1 < A: < a, awd iiwo sequences (q;„, n € N), 
{Pmn € N) of complex numbers such that: 

(a) the sequences (</3j(a„),n G N) and {ipk{l3n),n G N) are weZZ defined, 
(6) ft.i(a„) ^ anrf /i2(a„) 7^ 0, /or a// n G N, 

(c) /^n = </'j(an), /or all n G N, 

(d) lim„^oo (an,/3n, V'fc(/3n)) = (0,0,7). 

Proof. Proving the implication (m) (i) is easy. We now prove the implication [i) =^ [ii). 
By Lemma |4j there exists a sequence ((a„, 7„), n G N) in A'^ s.t. for all n G N we 
have: (1) /ii(a„) ^ 0, (2) /i2(a„) 7^ 0, (3) ri(a„,/3„) = 0, (4) r2(a„, /3„, 7„) = 0, (5) 
lim„_>.oo (ttn, /3n, 7n) = (0, 0, 7). FoUowing the proof of Lemma[7j we know that for n large 
enough the product (/?„ — (fi{an)) ■ ■ • (/?„ — <y3c(a„)) is zero. Therefore, from one of the 



8 



sequences (/?„ — 93i(a„)), . . . , (/3„ — ipc{an))i say the j-th, one can extract an (infinite) sub- 
sequence whose terms are all zero. Thus, without loss of generality, we assume that /3„ — 
(Pj{an) holds, for all n € N. Hence, for all n G N, we have /|(a„,7„) = r2(Q;„, /3„, 7„) = 0. 
Together with Equation ([3| and following the proof of Lemma [7j we deduce the desired 
result. □ 

Lemmas [6] and |8] immediately imply the following. 

Proposition 3. For all 7 e C, the following two conditions are equivalent. 

(i) (0,0,7) e lim(W^(ri,r2)) holds, 

(m) there exist integers j,k with 1 < j < c and 1 < k < a, such that ipj{0) — and 
V'fc(O) =7- 

Therefore, applying Puiseux theorem to ri and /I, then checking the constant terms of 
the series ipi, . . . ^ipb provides a way to compute all 7 £ C such that (0, 0, 7) is a limit point 
of W{ri, r2). Theorem [3] in Sections [4] states this principle formally for an arbitrary regular 
chain R in dimension one. 

Finally, one should also consider the case /ii(0) 7^ 0, /i2(0) = 0. In fact, it is easy to see 
that this latter case can be handled in a similar manner as the case /ii(0) = 0, /i2(0) = 0. 



4 Puiseux expansions of a regular chain 

In this section, we introduce the notion of Puiseux expansions of a regular chain, motivated 
by the work of [121 [T] on Puiseux expansions of space curves. 

Lemma 9. Let R =^ {ri, . . . ,rs_i} C C[Xi < • • • < Xs\ be a strongly normalized regular 
chain whose saturated ideal has dimension one. Recall that hji{Xi) denotes the product of 
the initials of polynomials in R. Let p > be small enough such that the set < \Xi\ < p 
does not contain any zeros of hn. Denote by Up :— {x — {xi, . . . ,Xs) G | < < 
p}. Denote by Vp{R) := vIr) nUp. Then we have W{R) D Up = Vp{K). Let R' := 
{primpart(ri), . . . , primpart(rs_i)}. Then Vp{R) = Vp{R'). 

Proof Let x G W{R)r]Up, then x e V{R) and x eUp hold, which implies that W{R)nUp C 
V{R) n Up. Let X e V{R) n Up. Since Up n V{hR) = 0, we have x S W{R). Thus 
F(i?) nUpC W{R) n Up. So W{R) nUp = Vp{R). similarly we have Vp{R) = Vp{R'). □ 

Notation 1. Let C C^ Denote limo(VF) {x = (xi, . . . , a;^) e | x e \im{W) and xi = 
0}. 

Lemma 10. Let R = {ri, . . . ,rs_i} C C[Xi < • • • < Xs]. Then we have \miQ{W{R)) = 
\ima{Vp{R)). 

Proof By Lemma[9) we have W(R) n Up{R) = Vp{R). On the other hand limo(W^(i?)) = 
limo(V(i?) n Up{R)). Thus \iYno{W{R)) = \\mo{Vp{R)) holds. □ 

Lemma 11. Let R be as in Lemma^ For 1 < i < s ~ 1, let di := deg(ri, X^+i). Then R 
generates a zero-dimensional ideal in C{{Xi))[X2, . . . , Xs]. Let V*{R) be the zero set of R 
in C{{X^))^~^ . Then V*{R) has exactly Yit^i points, counting multiplicities. 

Proof. It follows directly from the definition of regular chain, Bezout bound and the fact 
that C{{Xl)y^^ is an algebraically closed field. □ 
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Definition 1. We use the notions in Lemma 11 Each point in V*{R) is called a Puiseux 
expansion of R. 

Notation 2. Letrn^ \V*{R)\. WriteV*{R) = {$i,...,$™} mi/i = • • • , *r^(^i)); 

for i — 1, . . . , m. Let p > be small enough such that for 1 < i < m, I < j < s — I, each 
^iiXi) converges m < < p. We define V*{R) := U'^ilx e C | < |a;i| < p,Xj+i = 

Theorem 1. We have V*{R) = Vp{R). 

Proof. We prove tliis by induction on s. For i = 1, . . . , s — 1, recall that hi is the initial of 
r, . If s = 2, wc have 



(Xi,X2) -/^l(Xl)^(X2-$K^l))■ 



So V*iR) = Vp{R) clearly holds. 

Write i? = i?'U{r,_i}, X' = X2,...,X,_i, X = {X^, X' , X,)^x' = X2, . . . , 
X = {xi,x' ,Xs), and m' = |V"*(i?')|. For i = 1, . . . , to, let <i>j = where $■ stands 

for . . . , Assume the theorem holds for R' , that is Vp{R') = Vp{R'). For any 

i = 1, . . . , to', there exist ik G {1, . . . , to}, = 1, . . . , ds-i such that 

ds-l 

r,.i(Xi,X' = %,X,) = h,{X,) n (X, - <i>^;HXi)). (4) 

fe=i 

Note that V*{R) = U^'iUf^j/ {(X' = ^C^)}- Therefore, by induction hypothesis 

and Equation Q, we have 

v;{R) = yjti^1:=i {^\^^Up,x' = mx,),xs = ^;\x,)} 
= utY{^ I (^1,^') e v;{R'),Xs = K\xi)} 

= {x I e V7(i?'),r.-i(xi,a:',a;,) = 0} 

= {a; I e Vp{R'),rs^i{xi,x' ,x,) = 0} 

□ 

Theorem 2. LetV^g{R) {$ = (^i, . . . , € | ord($J) > 0,i = 1, . . . , s - 1}. 

T/ien we have 

limo(Ty(i?)) = U*ev|^(fl,){(Xi = 0,$(Xi = 0))}. 

Proof. By definition of V^>Q(i?), we immediately have 

limo(y:(i?)) - U*ey|„(fl){(^i = 0, - 0))}. 

Next, by TheoremjlJ we have V*{R) = Vp{R). Thus, we have limo(V';(i?)) = limo(Fp(i?)). 
Besides, with Lemma [lO| we have limo(W^(i?)) — \ima{Vp{R)) . Thus the theorem holds. □ 

Definition 2. LetV^^{R) := {$ = . . . , G | ord($J) > 0, j = l,...,s-l}. 

Let M = \V^q{R)\. For each = . . . , $,""^) e V^>o(i?); 1 < i < Af, we fcnow f/iai 

$1 G C((X]^)). Moreover, by Equation we fcnow that for j = 1, . . . , s— 1, is a Puiseux 

expansion ofrj{Xi, X2 — . . . , Xj — Xj^i). Let Q j &e t/ie ramification index of^l 



10 



and (T^'-J , Xj+i — ipj{T)), where tpj G C(r), be the corresponding Puiseux parametrization 
of ^j. Let Q be the least common multiple of {q.i, • • . , Let gf = ipl{T = T'"'/^*-'). 

We call the set <&r := {{Xi = T'-% X2 = g}{T), . .'. , X, = g^^-'^ (T)) , i = 1,...,M} a system 
of Puiseux parametrizations of R. 

Theorem 3. We have 

limo(VF(i?)) = ©fl(r = 0). 
Proof. It follows directly from Theorem [2] and Definition [2] □ 

5 Puiseux parametrization in finite accuracy 

In this section, we define the Puiseux parametrizations of a polynomial / G <C{X)\Y] in 
finite accuracy, see Definition [4] ^ 

For / e C(X)[y], we define the approximation f oi f for a given finite accuracy, see 
Definition 3 This approximation / of / is a polynomial in C[X,y]. In Section [7) we prove 
that in order to compute a Puiseux parametrizations of / of a given accuracy, it suffices to 
compute a Puiseux parametrization of / of some finite accuracy. 

In this section, we review and adapt the classical Newton-Puiseux algorithm to compute 
Puiseux parametrizations of a polynomial / G C[X, y] of a given accuracy. Since we do 
not need to compute the singular part of Puiseux parametrizations, the usual requirement 
discrim(/, 7^ is dropped. 

Definition 3. Let f = YT=o ^ For any t eN, let /(^) := J2^=o <^i^'- W'e call 

Z'"^' the polynomial part of f of accuracy r + 1. Now, let f — '^'^^Qai{X)Y^ e C{X)[Y]. 
For any t £ N, we call f^'^^ := X^^Lo <^i^^Y^ the approximation of / of accuracy t + 1. 

Definition 4. Let f e C{X)[Y], deg(/, F) > 0. Let a,T e N>o and g{T) = 'ZlZl^kT''. 
The pair [T" ^ g{T)) is called a Puiseux parametrization of / of accuracy r if there exists 
an irreducible Puiseux parametrization (T'' ,ip(T)) of f such that 
(i) a divides <;. 
(a) gcd{a,bo,...,br-i) = 1. 

(Hi) g{T''^'^) is the polynomial part of ip{T) of accuracy {<;/a)(T — 1) + 1. 

Note that if a — then g{T) is simply the polynomial part of ip{T) of accuracy t. 

We borrow the following notion from [14 in order to state an algorithm for computing 
Puiseux parametrizations. 

Definition 5. A C-tern^is defined as a triple t — {q,p,l3), where q and p are coprime 
integers, q > and j3 £ <C is non-zero. A C-expansion is a sequence tt = (^1,^2, • • •) of C- 
terms, where ti = {qi,pi, f3i), We say that tt is finite if there are only finitely many elements 
in TT. 

Definition 6. Let tt = (ti, . . . , ijv) be a finite C-expansion. We define a pair {T" ,g{T)) of 
polynomials in <C\T] in the following manner: 

• if N ^ 1, set a = I, g{T) = and 5n = 0, 

• otherwise, let a := Ilili?*.' := (P3llk=j+i1kj A < « < N), and 5^ 
Ci/gcd(a, ci, . . . , Civ) A < i < N). Set a :— a/gcd(a, ci, . . . , cjv) and g{T) :— 

^It is a simplified version of Duval's definition. 
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We call the pair {T" , g{T)) the corresponding Puiseux para-metrization of tt of accuracy 
6]y + 1. Denote by ConstructParametrization an algorithm to compute {T'^,g{T)) from tt. 

Definition 7. Let f e C{X)[Y] and write f as f{X,Y) J2Lo (E^o^^j^^) The 
Newton Polygon of f is defined as the lower part of the convex hull of the set of points 
in the plane such that ai,j ^ 0. 

Let / G C{X)[Y]. Next we present an algorithm, called Newton Polygon to compute the 
segments in the Newton Polygon of /. This algorithm is from R.J. Walker's book [2S\ . 
Newton Polygon (/, /) 

Input: A polynomial / e C{X)[Y]; a controlling flag /, whose value is 1 or 2. 
Output: The Newton Polygon of /. If / = 1, only segments with non-positive slopes 
are computed. If / = 2, only segments with negative slopes are computed. 
Description: 

— Write / as / = j:toUX)Y\ where b,{X) = Ej°1o«^j^'- 

— For < i < d, define 6i := ord(6i). 

— For < i < d, we plot the points Pi with coordinates («, Si); we omit Pi if Si = oo. 

— We join Pq to Pd with a convex polygonal arc each of whose vertices is a Pi and 
such that no Pi lies below the arc. 

— If / = 1, output all segments with non-positive slopes in the polygon; if / = 2, 
output all segments with negative slopes in the polygon. 

Next we present the specification of several other sub-algorithms which are necessary 
to present Algorithm [2] for computing Puiseux parametrization of some finite accuracy as 
defined in Definition |4l 

NewPolynomial(/,i,£) 

Input: /eC[X,r]; a C-term i = (g,p,/3); ^ G N. 
Output: A polynomial X"^ f{Xi,XP{P + Y)) e C[X,Y]. 

SegmentPoly(/. A) 

Input: / e C[X, F]; A is a segment of the Newton Polygon of /. 
Output: A quadruple (q,p,i,4>) such that the following holds 

— (7,p, ^ e N; (/) e C[Z]; q and p are coprime, q > 0. 

— For any (i, j) G A, we have qj +pi ~ £. 

— Let ^o min({^ | (^, j) G A}), we have ^ = E(»,j)eA a^^jZ^'-'o^'J . 
Theorem 4. Algorithm^ terminates and is correct. 

Proof. It directly follows from the proof of Newton-Puiseux algorithm in Walker's book [28] . 
the relation between C-expansion and Puiseux parametrization discussed in Duval's pa- 
per [131 , and Definitions [6] and [1] □ 
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Algorithm 1: NonzeroTerm(/, 7) 



Input: / e C[X, y]; / = 1 or 2 
Output: A finite set of pairs {t 
begin 
S:=0; 

for each A G Newton Polygon (/, /) do 
:= SegmentPoly(/, A); 
for each root ^ of (j) in <C do 

for each root j3 of If^ — ^ in C 
t := {q,P,l3); 
S := S[J{{t,i)} 



where t is a C-term, and ^ e N. 



do 



9 end 



Algorithm 2: NewtonPuiseux 



Input: / e C[X, Y]; a given accuracy r > e N. 

Output: All the Puiseux parametrizations of / of accuracy r. 

begin 

7r:=();5:={(7r,/)}; 
while 5 7^ do 

choose (tt*,/*) €S;S := 5 \ {(tt*, /* )}; 
if TT* = ( ) then I := 1 else / := 2; 
{T" ,g{T)) := ConstructParametrization(7r*); 
if deg{g{T),T) + 1 < r then 
C := NonzeroTerm (/*,/); 
if C = then 

output {T'',g{T)) //a finite Puiseux parametrization is 
found 

else 

for each {t = {p,q,/3),£) e C do 
TT** := TT* U (i); 

/** := NewPolynomial(/*,i,^); 

S := S'U{(7r**,/**)} 



16 

17 end 



output (T<^,5(T)) 
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6 Computing in finite accuracy 

Let / e C(J'(')[y]. In this section, we consider the foUowing problems. 

(a) Is it possible to use an approximation of / of some finite accuracy m in order to 
compute a Puiseux parametrization of / of some finite accuracy t? 

(6) If yes, how to deduce m from / and t? 

(c) Provide a bound on ni. 
Theorem [5] provides the answers to (a) and (5) while Lemma 15 answers (c). 

Lemma 12 ([15j. Let X = Xi, . . . , and Y ^ Yi, . . . ,Y^. Forgi,...,gs E C[[Y]], with 
OTd{gi) > 1, there is a C-algebra homomorphism (called the substitution homomorphism j 



c[[x]] c[[y]] 



f ^ f{9i{Y),...,gs{Y)). 

Moreover, if gi, . . . , g^ are convergent power series, then we have ^g{C{X)) Q C{Y_) holds. 

Definition 8 ([15]). Let f = ^a^^X^F'^ e C[[X,Y]]. The carrier of f is defined as 

carr(/) = {(M,i^) eN^ | a^.^O}. 

Lemma 13. Let f e C{X)[Y]. Let d := deg(/, F) > 0. Let q e N>o, p,£ €N and assume 
that q and p are coprime. Let /3 7^ G C. Assume that q,p,i define a line L : qj + pi = £ 
in (i, j) plane such that 

(a) There are at least two points (ji, ii) £ carr(/) and (j2, ^2) G carr(/) on L with ii ^ ii. 

(6) For any {j,i) € carr(/), we have qj +pi> £. 
Let fi := X^^ f(Xf, Xf{/3 + Yi)). Then, we have the following results 

{i) We have fi e C(Xi)[Yi]. 

(ii) For any given mi e N, there exists a finite number m € N such that the approximation 

of fi of accuracy mi can be computed from the approximation of f of accuracy m. 
(Hi) Moreover, it suffices to take m — [^ "i+^ J . 

Proof. Since q> {) holds, we know that oid{X'l) = q > holds. We also have fiXf, Xf{f3 + 
Yi)) e C(Xi)[ri]. Let f{X,Y) ^to (E^to^^J^') Y\ Then we have Yi) = 

Ef=o (E^o^^j^i'^''^^*"^^) (-^ + ^i)'- Since for any {j,i) e carr(/), we have qj + pi > e, 
the power of Xi cannot be negative. By Lemma 12 we have /i e C(Xi)[Yi]. That is (i) 
holds. 

We prove (ii). We have 



/i(Xi,Yi) mod (Xf^) 



Since q € N>o and mi, £ and i are all finite, we know that j has to be finite. In other words, 
there exists a finite m such that the approximation of /i of accuracy toi can be computed 
from the approximation of / of accuracy to. That is, (ii) holds. 

Since the first toi terms of /i depends on the j-th terms of /, which satisfies the 
constraint qj+pi~£< mi , we have j < (™i+^^)~p' < (™i+^) ^ j^^^ ^-^^ maximum 

of these j's . Now we have m' — 1 < ("'1^+^) ^ Since m' is an integer, we have to' < [ ("^^^^^ J 

holds. Let m ^ [ilT^J . Next we show shat toi > 1 implies that to > 1 holds. If there is 
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at least one point G L such that j > 1, then we have £ > q, which imphes to > 1. If 
the j-coordinates of all points on L is 0, then q = I and £ = 0, which implies also m > 1. 
Thus (iii) is proved. □ 

Remark 1. We use the same notations as in the previous Lemma. In particular, let 
f{X,Y) ^to {ET^oa^jX^) <^^d h ■= X^' f{Xf, XfijS + Y,)). For a fixed term 
aijX^Y'' of f, it appears in fi as 



a.,,Xr^P'-'{P + Y,y ^}_^[[ l]P^-'a.,,Xr^P''-' ) Y,'. 

fc=0 



For two fixed terms aij-^X^^Y^ and Oi^j^X^'^Y'^ of f with ji < j2, since qji + pi — £ < 
132 + pi — £, we know that for any fixed k, Oi^j^X^'^Y^ always contributes strictly higher 
order of powers of Xi than Oij-^X^^Y^ in f\. 

Remark 2. Let f{X,Y) := (E^o«^j^^) ■ For < i < d, let a^.J> he the first 

nonzero coefficient among {a^ j|0 < j < 00}. We observe that the Newton polygon of f is 
completely determined by ai,j- , <i < d. 

Theorem 5. Let f e C{X)[Y]. Let t e N>o. Let a e N>o and g{T) T.lZlbkT''. 
Assume that (T'^ , g{T)) is a Puiseux parametrization of f of accuracy r. Then one can 
compute a finite number m S N such that (T"' , g{T)) is a Puiseux parametrization of ac- 
curacy T of the approximation of f of accuracy m. We denote by AccuracyEstimate an 
algorithm to compute such m from f and r. 

Proof. Let /o := /, ^0 ^ and Yq Y. For i = 1,2,..., Newton-Puiseux's algorithm 
computes numbers qi,pi,£i, Pi and the transformation 

/,:=xr'7.-i(^f,xr(A + y.)) 



such that the assumption of Lemma 13 is satisfied. 

By Lemma |13[ we know that for any i, a given number of terms of the coefficients of 
fi in Yi can be computed from a finite number of terms of the coefficients of /i„i in Yi^i. 
Thus for any i, a given number of terms of the coefficients of fi in Yi can be computed from 
a finite number of terms of the coefficients of / in Y. 

On the other hand, the construction of Newton-Puiseux's algorithm and Remark [5] tell 
us that there exists a finite M, such that a and all the terms of ^(T) can be computed from 
a finite number of terms of the coefficients of fi in Yi, i — 1, ... , M . 

Thus we conclude that there exists a finite number m e N such that {T'^,g{T)) is a 
Puiseux parametrization of accuracy r of the approximation of / of accuracy to. 

Next we show that there is an algorithm to compute to. We initially set to' := r. Let 

/o := J2i=o (SjLo^jj^"') That is, /o is the approximation of / of accuracy to' + 1. 

We run Newton-Puiseux's algorithm to check whether the terms ak,m'X"^ Y^ , < k < d, 
make any contributions in constructing the Newton Polygons of all fi. If at least one of 
them make contributions, we increase the value of to' and restart the Newton-Puiseux's 
algorithm until none of the terms (Xf^ rn' 

X"^ Y^, < k < d, makes any contributions in 
constructing Newton Polygons of all fi. By Remark [l} we can set to := to'. □ 
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Lemma 14. Let d,T E N>o. Let Uij, 0<i<d, 0<j<T, and bk, < k < t be symbols. 
Write a = (ao^o, ■ ■ • ,ao,r-i, 

■ . ■ ,ad.o, ■ ■ ■ ,ad,T-i) and b = (60, • ■ • , &r-i)- Let f{a,X,Y) = J2i=a (EJ=o Oij"^'') ^ 
C[a][X,r]. Letg(b,X)-El;o 

bfeX'^' e C[b][X]. Let p /(a,X,r = .g(b,X)). Let Ffc := cocff(p, X*^), < fc < t - 1}, 
and F {i^o, • • • , ^r-i}- 7"/ien under the order a < b and bg < bi < ■ ■ ■ < br-i, F forms 
a zero- dimensional regular chain in C(a)[b] with main variables [bg, 61, ... , and main 

degrees (d, 1, . . . , 1). In addition, we have 

• Fo = J2t=o «i.o^o ^"■'^ 

• init(Fi) = • • • = init(F.r-i) = « ' 04,0^0"^ ■ 

Proof. Write p — X^iLo {j2^=o o-i.jX^^ (Sfe=o ^kX^^ as a univariate polynomial in X. 

Observe that Fq = X)iLo "^^o^o- Therefore Fq is irreducible in C(a)[b]. Moreover, we have 
mvar(_Fo) = ^0 a-nd mdeg(i^o) = d. 

Since d > 0, we know that ai^ (^X]fe=o ^kX^^ appears in p. Thus, for < /c < r, bk 
appears in Fk- Moreover, for any k > 1 and i < k, bk can not appear in Fi since bk and 
X'' are always raised to the same power. For the same reason, for any i > 1, 6^ cannot 
appear in Fk, for 1 < fc < r. Thus {Fq, . . . ,Ft-^i} is a triangular set with main variables 
(60, bi, . . . , 6r-i) and main degrees (d, 1, . . . , 1). 

Moreover, we have init(_F'i) = ■ ■ ■ — init(i^T-i) = X]i=i * 'Qi,o&o^^i which is coprime with 
Fq. Thus F = {Fq, . . . , Fr-i} is a regular chain. □ 

Lemma 15. Let f = J2Lo {'ET=o "-^.j^^) ^' Assume that deg(/, F) > 

and f is general in Y. Let (p{X) = X^fcLo ^^"'^'^ ^ '^[[^1] such that f{X,(p{X)) = holds. 
Let r > e N. Then "generically" , bi, < i < t, can be completely determined by 
{a,j I < i < d,0 < j < r}. 

Proof By f{X, Y) ^ 0, we know that f{X, Y) ^ mod (X^). Therefore, we have 

i=0 \j<T J \k<T ) 

Let p = Eto i^,<r^^■J^') {EkKrbkX'^y. Let := {coeS{p,X^), < i < r}, and 
F := {i^o, ■ ■ • , -FV-i}- Since / is general in Y and /(X, = 0, there exists i* > 



such that Oi* ,0 7^ 0. By Lemma 14 we have Fq = X)f=o '^^ o^o- Thus 69 can be completely 



determined by ai^, < i < d. In order to completely determine 61, ... , 67—1, it is enough to 
gurantee res(i^o, Fi, 69) 7^ holds. Therefore the values of 6^, < /c < r can be completely 
determined from almost all the values of Oij, 0<i<d, 0<j<T. □ 



7 Accuracy estimates 

Let R :— {ri{Xi,X2), . . . ,rs~i{Xi, . . . ,Xs)} C C[Xi < • • • < Xg] be a strongly normalized 
regular chain. In this section, we show that to compute the limit points of W{R), it suffices 
to compute the Puiseux parametrizations of R of some accuracy. Moreover, we provide 
accuracy estimates in Theorem [6] 
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Lemma 16. Let f = adiX)^^ H h ao{X) e C{X)[Y], where d > and ad{X) ^ 0. 

For < i < d, let 5i :— ord(ai). Let k :— min((5o, . . . ,Sd)- Let f f / X^ . Then we have 
f £ C{X)[Y] and f is general in Y. This process of producing f from f is called "making 
/ general" and denote by MakeGeneral an operation which produces f from f. 

Proof. Since k = niin(5o, . . . , Sd), there exists i, I < i < d, such that k = Si. Moreover, for 
all 1 < J < rf, we have Sj > k. Thus for every such i, we have OTd{ai{X) / X'') = and 
aj{X)/X'' e C{X), 0<j<d. This shows that / G C{X)[Y] and / is general in Y. □ 

The following lemma shows that computing limit points reduces to making a polynomial 
/ general. 

Lemma 17. Let f G C{X)[Y], where deg(/, y) > 0. Assume that f is general in Y. Let 
p > be small enough such that f converges in \X\ < p. Let Vp{f) :— {{x,y) S | < 
\x\ < p, fix, y) = 0}. Then we have limo(y,(/)) = {(0, G | /(O, y) = 0}. 

Proof. Let {X — T''',Y — Lpi{T)), 1 < i < c < d, he the Puiseux parametrizations of 
/. By Lemma |9] and Theorem [sj we have limo(V^p(/)) = U^^i{(0,?;) G | y = (y3,(0)}. 
Let (X = T'^' , gi(T)), i = 1, . . . ,c, be the corresponding Puiseux parametrizations of / of 
accuracy 1. By Theorem [sj there exists an approximation / of / of some finite accuracy 
such that (X = T"^' , gi{T)), i — 1, . . . ,c, are also Puiseux parametrizations of / of accuracy 
1. Thus, we have ipi{0) = gi{0), i — 1, . . . ,c. Since / is also general in Y, by Theorem 2.3 
of Walker [2H], we have U,=^i{(0,?/) e \ y = g^iO)} = {(0,?/) G | /(0,y) = 0}. Since 
/(O, y) = /(O, y), the Lemma holds. □ 

Lemma 18. Let a{Xi, . . . , Xs) E C[Xi, . . . , X^]. Let gi = J^T^o ^^J^^ ^ '^(^) ■ ^"^^ 
a{gi, . . . ,gs) as X]fe°=o ^kP*' . To compute a given bk, one only needs the set of coefficients 
{c.j \ l<i<s,0<j<k}. 

Proof. We observe that any Ci^, where j > fc, does not make any contribution to 6^.. □ 

Lemma 19. Let f = ad{X)Y'^ -\ \-ao{X) G C{X)[Y], where d>0, and ad{X) ^0. Let 

S := OTd{ad{X)). Then "generically", a Puiseux parametrization of f of accuracy r can be 
computed from an approximation of f of accuracy t -\- 5. 

Proof. Let / := MakeGeneral(/). Observe that / and / have the same system of Puiseux 



parametrizations. Then the conclusion follows from Lemma [16^ and Lemma 15 □ 



Theorem 6. Let R := {ri(Xi, X2), . . . , rs_i(Xi, . . . , X,)} C C[Xi < ■■■ < X^] be a 
strongly normalized regular chain. For \ <i < s — 1, let hi init(ri), di :— deg{ri, Xi^i) 
and Si :— ord(/ii). We define fi, 2 < i < s — 1, and <;_,■, Tj, (pj(Tj), 1 < j < s — 2, as follows 

• Let (Xi — T^'^ , X2 — ipi(Ti)) be a Puiseux parametrization of ri(Xi, X2). 

• Let f, := r,(Xi = T^\ X2 - MTi), .■.,X,^ (^^-i^^i), 
Xi+i). 

• Let {Ti^i = Tl'-.Xi^i — ifi(Ti)) be a Puiseux parametrization of fi. 



Then we have the following results: 
(i) Let To := Xi, for < i < s ~ 2, define gi{Ts^2) '■= yj^|='+i ^ then we have Ti 

gi{Ts-2)- 
[li) We have G C(r,_2)[^.]. 



17 



(in) There exist numbers ti, . . . , Ts_2 G N such that in order to make fs-i general in Xs, it 
suffices to compute the polynomial parts of ipi of accuracy Tf , 1 < i < s — 2. Moreover, 
if we write the algorithm AccuracyEstimate for short as 9, the accuracies can he 
computed in the following manner 

- let Ts-2 (Ilfc^i ^k)Ss-i + 1 

- let n-i := ma.x{6{fi,n), (JXk'Ji'ik)Ss-i + I), for s - 2 > i > 2. 

(iv) Generically, for 1 < i < s — 3, we can choose Ti = (0^=1 '>fe)(X]fc=2 ^i) + 1- 
(w) The indices (^k can be replaced with dk, k = 1, . . . , s ~ 2. 

Proof. We prove (i) by induction. Clearly [i) holds for i = s ~ 2. Suppose it holds for i. 
Then we have 

Therefore (i) holds also for i — 1. So (i) holds for all < i < s — 2. 
Note that 

fs-i - r,^iiXi^T',\X2^ipi{Ti),-.-, 

= rs-i{Xi=ga{Ts-2),X2^ipi{gi{Ts-2)),---, (5) 

Xs-2 = iPs-3igs-3iTs~2)), 
Xs-1 = ips-2igs-2iTs-2)), Xs) 

Since 0Td{gi{Ts^2)) > for all < « < s — 2, by Lemma [l2| (m) holds. 
Note that go{Ts-2) = tI!^^'^''- Since ord(/is-i(^i)) = Ss^i, we have 

/s-2 \ 



ord(/i,_i(Xi = 5o(T,_2))) =11'^^ '^^-1- 



Let Ts-2 '■= {Y[k=^i^k)Ss-i + 1- By Lemma 16 to make fs-i general in Xg, it sufBces to 



compute the polynomial parts of the coefficients of fs-i of accuracy Ts-2- 

By Lemma jTsj and Equation ([s]), we need to compute the polynomial parts of (pi{gi{Ts-2)), 
1 < z < s — 2, of accuracy Ts-2- Since ord((?i(Ts_2)) = Ilfe^i+i'^fc' achieve this accu- 
racy, it's enough to compute the polynomial parts of ipi of accuracy (111=1 'ik)Ss-i + 1, for 
1 < i < s - 2. 

On the other hand, since /, = r,(Xi = ,^2 = MTi), ■ • ■ = (/7,_i(r,_i), X,+i) 
and (Ti_i = T^^Xi^i = tfi{Ti)) is a Puiseux parametrization of fi, by Theorem [5] and 
Lemma |l8[ to compute the polynomial part of (pi of accuracy r^, we need the polynomial 
part of ipi-i of accuracy 0{fi,Ti). 

Thus, take Ts-2 := (nfe=i ^fe)'^s-i + 1 and lp^-l = max(6'(/i, n), (01=1 ?fc)^s-i + 1) for 



2 < i < s — 2 will guarantee fs-i can be made general in Xs- So (Hi) holds. By Lemma 19 
generically we can choose 6{fi, Ti) = Ti + (111=1 '^k)^i, 2 < i < s — 2. Therefore (iv) holds 
Since we have ^ dk, 1 < k < s — 2, {iv) holds. □ 



8 Algorithm 

In this section, we provide a complete algorithm for computing the non-trivial limit points 
of the quasi-component of a one-dimensional strongly normalized regular chain based on 
the results of the previous sections. 
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Algorithm 3: LimitPointsAtZero 



Input: R := {ri{Xi,X2), . . . ,rs_i(Xi, . . .,Xs)} C C[Xi < ■ ■ ■ < X,], s > 1, is a 

strongly normalized regular chain. 
Output: The non-trivial limit points of W{R) whose -coordinates are 0. 
1 begin 

let S := {(To)}; 

compute the accuracy estimates ti, . . . ,Ts_2 by Theorem[6j let Ts„i = 1; 
for i from 1 to s — 1 do 
y:=0; 
for $ e S* do 

f, := r^{Xi = $1, . . . ,Xj = 
if i > 1 then 
I let 6 := ord(/„T,_i); let /, := /./Tf_i; 

E Newton Pu iseux(/i, Tj); 
for (r,_i = ^{T,),X,+, = (^(T,)) e S do 
I S' := 5' U {ci>(T,_i = U (^(T,))} 

5* := S" 

if 5 = then return else 
|_ return eval{S,Ts-i = 0) 



16 end 



Algorithm 4: LimitPoints 



Input: A strongly normalized regular chain 
i?:={ri(Xi,X2),...,r,_i(Xi,...,X,)}cC[Xi < 
Output: All the non-trivial limit points of W{R). 
1 begin 



< A,], s > 1. 



let hii := init(i?); let L be the set of roots of Hr; 

5:=0; 

for a G L do 

Ra := R[Xi = Xi + a)- 
Sa ■= LimitPointsAtZero(i?a); 

update Sa by replacing the first coordinate of every point in it by a; 
S:^SUSa 



return S 



10 end 
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Remark 3. Note that line 9 of Algorithm\^ computes Puiseux parametrizations of fi of 
accuracy t^. Thus {(j){Tj) , Lp{Tj)) at line 10 cannot have negtive orders. 

If the D5 principle is applied to Algorithms^ and [7} the limit points ofW{R) can be 
represented by a finite family of regular chains. 

Proposition 4. Algorithm^is correct and terminates. 

Proof. It follows from Theorem [3] Theorem [5] Theorem [6] and Lemma |17| □ 

9 Experimentation 

We have implemented Algorithm [4] of Section |8j which computes the limit points of the 
quasi-component of a one-dimensional strongly normalized regular chain. The implementa- 
tion is based on the library RegularChains and the command algcurves[puiseux] of Maple. 
The code is available at http : //www. orcca. on. ca/~cchen/ACM13/LimitPoints .mpl. This 
preliminary implementation relies on algebraic factorization, whereas, as suggested in |14j . 
applying the D5 principle, in the spirit of triangular decomposition algorithms, for in- 
stance [5], would be sufficient when computations need to split into different cases. This 
would certainly improve performance greatly and this enhancement is work in progress. 

As pointed out in the introduction, the computation of the limit points of the quasi- 
component of a regular chain can be applied to removing redundant components in a Kalk- 
brener triangular decomposition. In Table [TJ we report on experimental results of this 
application. 

The polynomial systems listed in this table are one- dimensional polynomial systems 
selected from the literature |S1[S]. For each system, we first call the Triangularize command 
of the library RegularChains, with the option "'normalized='strongly', 'radicar = 'yes"' . For 
the input system, this process computes a Kalkbrener triangular decomposition TZ where the 
regular chains are strongly normalized and their saturated ideals are radical. Next, for each 
one-dimensional regular chain R in the output, we compute the limit points lim(Vl^(i?)), thus 
deducing a set of regular chains i?i , . . . , i?e such the union of their quasi-components equals 
the Zariski closure W{R). The algorithm Difference [S is then called to test whether or not 
there exists a pair R, R' of regular chains of TZ such that the inclusion W{R) C W{R') 
holds. 

In Table[T] the column T and #(T) denote respectively the timings spent by Triangularize 
and the number of regular chains returned by this command; the column d-1 and d-0 
denote respectively the number of 1-dimensional and 0-dimensional regular chains, whose 
sum is exactly #(T); the column R and 7^(R) denote respectively the timings spent on 
removing redundant components in the output of Triangularize and the number of regular 
chains in the output irredundant decomposition. As we can sec in the table, most of the 
decompositions are checked to be irredundant, which we could not do before this work 
by means of triangular decomposition algorithms. In addition, the three redundant 0- 
dimensional components in the Kalkbrener triangular decomposition of system f-744 are 
successfully removed. Therefore, we have verified experimentally the benefits provided by 
the algorithms presented in this paper. 
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Table 1: Removing redundant components. 



Sys 


T 


#(T) 


d-1 


d-0 


R 




f-744 


14.360 


4 


1 


3 


432.567 


1 


Liu-Lorenz 


0.412 


3 


3 





216.125 


3 


MontesSS 


0.072 


2 


2 





0.064 


2 


Neural 


0.296 


5 


5 





1.660 


5 


Solotareff-4a 


0.632 


7 


7 





32.362 


7 


Vermeer 


1.172 


2 


2 





75.332 


2 


Wang-1991c 


3.084 


13 


13 





6.280 


13 



10 Concluding remarks 

We conclude with a few remarks about special cases and a generalization of the algorithms 
presented in this paper. 

Reduction to strongly normalized chains. Using the hypotheses of Lemma [Sj we observe 
that one can reduce the computation of lim(VF(i?)) to that oi lim{W (N)) . Indeed, under the 
assumption that sat(i?) has dimension one, both lim(Vl^(_R)) and \im{W{N)) are finite. Once 
the set lim(W^(7V)) is computed, one can easily check which points in \ini{W{N)) do not 
belong to W{R) and then deduce \im{W{R)). This reduction to strongly normalized regular 
chains has the advantage that /ijv is a univariate polynomial in C[Xi], which simplifies 
the presentation of the basic ideas of our algorithms, see Section |3] However, it has two 
drawbacks. First the coefficients of N are generally much larger than those of R. Secondly, 
lim(W^(A^)) may also be much larger than lim(M^(i?)). A detailed presentation of a direct 
computation of lim{W{R)), without reducing to lim(Vl^(7V)), will be done in a future paper. 

Shape lemma case. Here, by reference to the paper [2] (which deals with polynomial 
ideals of dimension zero) we assume that, for 2 < i < e, the polynomial involves only 
the variables Xi,X2,Xi and that deg{ri, Xi) — 1 holds. In this case, computing Puiseux 
series expansions is required only for the polynomial of R of lower rank, namely ri . In this 
case, the algorithms presented in this paper are much simplified. However, for the specific 
purpose of solving polynomial systems via triangular decompositions, reducing to this Shape 
lemma case, via a random change of coordinates, has a negative impact on performance and 
software design, for many problems of practical interest. In contrast, the point of view of 
the work initiated in this paper is two-fold: first, deliver algorithms that do not require any 
genericity assumptions; second develop criteria that take advantage of specific properties of 
the input systems in order to speedup computations. Yet, in our implementation, several 
tricks are used to avoid unnecessary Puiseux series expansions, such as applying the theorem 
(see |15| p. 113) on the continuity of the roots of a parametric polynomial. 

Handling the case where sat(7?) has dimension greater than 1. From now on, sat(i?) 
has dimension s — e > 2. We use the notations of Lemma [5] and recall that each point of 
lim(M^(i?' U Te)) is in particular a point of lim{W{R')). Since we know how to compute 
\im{W{R')) when R' consists of a single polynomial, we assume, by induction that a trian- 
gular decomposition of \im{W{R')) has been computed in the form W{Ri) U • • • U W{Rf) 
for regular chains i?i , . . . , i?/ . 

We observe that a point p e lim(M^(i?')) can be "extended" to apoint of lim(Vl^(i?' U rg)) 
in two ways. First, if p does not cancel the initial of (which can be tested algorithmically), 
then, by applying the theorem (see again [15] p. 113) on the continuity of the roots to r^, 
we extend p with the Xg-roots of r^, after specializing {Xi, . . . ,Xs-i) to p. From now 
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on, we assume that p cancels the initial he of r^. In this case, we compute a truncated 
Puiseux parametrization about p using the regular chain Ri such that p £ W{Ri) holds. 
After substitution into the polynomial r^, we apply Puiseux theorem and compute the limit 
points of W{R' U Te) extending p, in the manner of the algorithms of Section [sj 

There are new challenges, however, w.r.t. to the one-dimensional case. First, parametriza- 
tions may involve now more than one parameter. When this happens, one should use 
Jung-Abhyankar theorem [25] instead of the Puiseux theorem. The second difficulty is that 
lim(M^(i?')) n V{he) may be infinite. This will not happen, however, if sat(i?') has dimen- 
sion at most 2 and is regular w.r.t. sat(i?'). This second assumption can be regarded 
as a genericity assumption. Thus the algorithms presented can easily be extended to di- 
mension two, under that assumption, which can be tested algorithmically. Overcoming in 
higher dimensions this cardinality issue with lim(W^(i?')) H V{he), requires to understand 
which "configurations" are essentially the same. Since lim(VF(i?)), as an algebraic set, can 
be described by finitely many regular chains, this is, indeed, possible and work in progress. 
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